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bustion  process  in  a  turbulent  flow; 
collection  of  articles].  Moskva,  Ob- 
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The  calculation  of  the  proc*?8ses  taking  place  In  i*ajnjet  com¬ 
bustors  is  a  complex  problem. 

Both  here  ([5]  and  [?])  and  abroad  ([!]#  [12],  and  [13]) 
attempts  have  been  made  to  develop  such  a  calculation  method.  . 

Usually,  the  poaltioo  of  the  flame,  the  pressure  drop,  and  the 
change  in  completeness  of  combustion  along  a  combustor  with  a  flame 
propagating  from  a  point- ignition  source  are  determined.  In  a  num¬ 
ber  of  works  (e.g.,  [7]),  this  problem  is  worked  out  for  a  flame 
front  of  finite  width. 

The  most  precise  mathematical  solution  has  been  obtained  by 
A.  V.  Talantov,  who  worked  out  the  problem  by  Introducing,  as  did 
other  authors,  a  number  of  sln^jllfications  (a  homogeneous  mixture, 
point- ignition  source,  constant  parameters  of  turbulence  across  and 
along  the  combustor,  etc.}.  However,  these  siBq;>lifications  do  not 
make  it  possible  completely  to  overcome  the  following  difficulties. 

(1)  Closing  the  system  of  basic  equations,  which  is  neces¬ 
sary  to  obtain  the  solution,  requires  the  application  of  axjxlllary 
conditions.  These  are  not  always  chosen  on  sufficiently  substanti¬ 
ated  grounds. 

(2)  The  precise  dependence  of  the  velocity  of  turbulent  flame 
propagation  on  the  initial  parameters  (u*,  Ug#  etc.)  is  not  known. 

(3)  There  is  no  uniform  view  about  the  way  in  which  the  width 
of  the  combustion  zone  of  a  turbulent  flame  (or  comibustlon  time)  should 
be  determined.  Experimental  data  on  the  determination  of  the  width 

of  the  combustion  zone  are  few  and  contradictory. 

In  :hls  connection,  it  is  felt  that  the  correct  solution  of 
the  problem  must  be  found.  This  is  the  purpose  of  the  present  work. 


1.  Survey  of  the  Literature 


The  solution  of  the  problem  of  flame  propagation  in  a  cylin¬ 
drical  tube  with  ignition  from  a  point  source  was  attempted  by  Ya.  B. 
Zel‘dovlch,  Q.  Tzyan  [H.  Tslen?],  G.  I.  Taganov,  A.  C.  Scurlock, 
fT.  Fabry,  etc.  These  authors  attempted  to  solve  the  problem  for  a 
^.aminar  flame  (the  thickness  of  flame  front  being  infinitesimal  and 
the  velocity  of  flame  propagation  low  in  ccmiparison  with  that  of  the 
incident  flow).  In  other  words,  the  flame  was  considered  as  a  zone 
of  combustion  products  separated  from  the  fresh  mixture  by  an  infi¬ 
nitely  thin  surface  (the  flame  front)  at  which  heat  is  momentarily 
emitted  as  a  result  of  the  combustion  reaction;  the  temperature  then 
rises  abruptly  by  several  times. 
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Such  a  simplification  is  too  \xnreal.  In  practice  the  width  of 
the  ccaahustion  zone  may  be  substantial  and  isast  not  be  disregarded 
(see  [4]  and  [8]). 

Such  a  problem  was  worked  out  for  a  eosdjustion  zone  of  finite 
width  by  B.  P.  Skotnikov  and^  later,  by  A.  Y,  TalantoT  [73. 


Fig.  1.  Diagram  of  a  flame  in  a  tube 


In  calculating  a  combustor  the  problem  is  usually  stated  as 
follows.  The  combustor  has  a  constant  cross  section,  xne  mixture 
is  Ignited  at  a  point  on  the  axis  of  the  flow.  Before  eateriiig  the 
coMibustor,  the  flow  is  one- dimensional.  The  static  pressure  at  the 
combustor’s  cross  section  is  assumed  to  be  constant  (Pig.  l).  The 
laws  gover  ilng  changes  of  temperature  and  velocity  in  the  ccmibustiow 
zone  in  the  transverse  direction  are  considered  to  be  given.  It  is 
necessarj"  to  find  a  method  for  calculating  all  the  parameters  of  the 
flow  in  an  arbitrary  cross  section  /.  The  following  values  are  de¬ 
termined;  the  velocities  of  the  fresh  mixture  (wc)  and  of  the  com¬ 
bustion  products  (wn),  the  temoeratures  of  the  fresh  mixture  (Tc)  axid 
of  the  combustion  products  (Tjil,  pressure  (px)^  the  ordinate  of  the 
Initial  flame  front  (yc)#  and  the  ordinate  of  the  end  of  combustion 
zone  (yn).  The  equations  given  below  are  used  for  the  solution*. 


The  energy  equation  for  a  fresh- mixture  flow  is 


0  T  ^ 


4' 


Where  Tg  ajyi  Wg  are  the  tesperature  and  flow  velocity  at  the  combust ox 


tFe  basic  equations  absolutely  necessary  for  any  calculation  of 
a  combustor  are  given  below. 
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entrance j,  respectively;  or,  after  transformation,  it  is 


«h«r« 


’e' 


^  is  the  Mach  number  at  combustor  entrance. 

The  equation  for  the  adiabatic  process  is 


7*  V/>o/ 


or 


wher* 


k~\ 


(2) 


is. 

P9 


The  energy  equation  for  the  flow  of  c<M!nbustion  products  is 


B 


or 


M  k  =i+_l- 

/q  Wo  9 


(3) 


The  equation  of  conservation  of  mass  for  the  entire  flow  is 


(/"o  ~  ^«)  +  P.w.i, -f.  J  dF, 


(4) 


where  P^,  and  Pn  are  the  densities  of  the  fresh  mixture  and  the 
'jombustion  products,  respectively  [sic],  in  the  initial  cross  section; 
Hg  and  Wg  are  the  density  and  flow  velocity  of  the  mixture  in  the 
combustion  zone;  and  Pn  are  the  cross-section  areas  of  the  com¬ 
bustor  and  of  the  flow  of  combustion  products;  and  Po  «  P^  -  (Fn  +  Pg) 
represents  the  area  of  the  entire  flow  minus  the  flow  of  fresh  ffilxtu2?e; 


«  3  . 


or 


vher® 


1=  *  «c(l  "  -Vn  +  M’ 

te  "a  J 


/  ■■  /  ^"-  •  It  --  '**'  •  T  3= 

/«  c  '  /n  -  ..  •  '^s—  ...  «  y 

,  '  0  /^O  « 0  /^6 


{•'^) 


The  equation  of  quantity  of  motion  is 


f 

PJ’,  +  '-.rf.  +  ,,  f.)  i-  ,tF. 


(6) 


or,  after  making  simple  transformations. 


MTwro 


i 

«> 


+  1- 


Z  .  /’I  /»«. 
foti’d 
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(7) 


The  equation  of  conservation  of  mass  of  the  fresh  mixture  is 


f'c(.  Wc{,.l,"Pc/{/^0  (8) 


Where  Pqq  is  the  average  density  in  a  portion  between  neighboring 
cross  s^tions,  and  AS  is  the  area  of  the  flame  front  between  those 
sections,  or,  in  a  dimensionless  form. 


0  an.^  0 

h  •• 
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In  order  to  solve  this  systexa  of  equations  It  is  necessary  to 
introduce  additional  conditions.  In  works  [7]  and  [11]  these  condi¬ 
tions  are  not  similar  and,  therefore,  the  final  resists  differ. 

In  A,  V.  Talantov‘s  work  the  eqxiations  needed  are  found  from 
the  condition  of  "equality  of  the  requia?ed  and  available  time  for 
the  streamtube,  "  IRie  available  time,  when  the  asotion  is  along  the 
line  of  flow  00*  (see  Fig.  l),  is 


The  tiB^  required  to  complete  the  combustion  process  can  be  found 
when  the  law  of  burnout  of  mixture  in  the  streamtube  is  known  in 
terms  of  time: 


(to) 

where  Ax  is  the  distance  between  neighboring  cross  sections,  is 
the  average  speed  of  motion  of  the  mlxtui?©  between  those  sections, 
tn  is  the  total  time  of  combustion  ©f  the  mixttsre  in  the  strsamtube, 
and  t3(z-l)  is  the  time  during  which  the  mixture  remains  in  the  zone, 
beginning  with  the  moment  it  crosses  the  initial  boundary  of  the  flame. 

In  addition,  the  author  makes  use  of  an  auxiliary  approximate 
condition  first  introduced  by  6.  Tzyan  [9].  Equations  (l)  and  (3) 
are  en^loyed: 


^  la—  j  I) 

or 

Cli) 

Finally,  the  portion  Qb  of  the  heat  libf rated  at  any  cross 
section  as  related  to  the  total  heat  available  la  the  mixture 
'burnout)  can  be  found  from  the  relation 


(12) 

Qo  Q« 

where  and  L  are  the  total  heat  content  in  the  mixture  at  cross 
section  i  and  at  the  initial  cross  section,  respectively. 
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one  obtains 


After  transformation. 


+'~l<y.(X,~l)  +  >: 


/<  /-  ^ 

,  (13) 

/»  /a  ^  i  ’ 


The  system  is  closed  when  the  laws  of  change  of  the  param¬ 
eters  of  the  coBibustion  zone  are  toiown  and  the  integrals  in  the 
equations  (5),  (7),  and  (13)  aa§  evaluated,  Talahtov  assumes  that 
the  law  of  change  of  velocity  and  temperature  in  the  zone  is  linear. 

^e  most  questionable  point  in  Talantov's  calculation  method 
is  his  introduction  of  a  characteristic  corabustion  time  in  a  turbulent 
flow;  the  combustion  time  is  selected  on  the  basis  of  the  author’s 
own  data.  [8] 

The  accuracy  of  both  the  jaethods  and  the  results  of  measuring 
tgji  are  debatable.  This  matter  is  discussed  in  detail  below. 

Virtually  the  same  equations  C(l)^  (^)s  and  (8)3  are  con¬ 

sidered  by  B.  P.  Skotnikov, 

The  following  additional  conditions  were  used  to  solve  the 
equations; 

(1)  The  connection  between  Uq  and  uw®  (where  una  is  the  ve¬ 
locity  of  combustion  products  along  the  ccmbustlon  axis)  was  determined 
in  the  form  of  the  relation  ' 


where  t  is  the  degree  of  preheating. 


fa 

(2)  The  law  of  change  of  the  velocity  profile  along  the  cross 
section  of  combustor  (and  of  the  density  profile)  was  defined  as  para¬ 
bolic. 
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(3)  The  dependence  -  tig  ojt‘  the  width  of  the  coaibustlon 

zone  in  a  rMial  direction  was  taken  in  the  form  6^  «  f(t^c)> 

which  experimental  data  were  utilized. 

Here  i=  the  radius  of  the  Initial  boundary  of  the  flame; 
tig  s  the  radius  of  the  boundary  of  the  flame  eosibustion 
products, 

Q^e  calculation  by  P.  Skotnikov  is  also  not  entirely  cor¬ 
rect,  The  definitions^  «  f(tQ)  is.  Incorrect,  since  depends  not 
on  tig  but  on  S,  which  is  the  distance  from  the  point  of  Ignition  of 
the  fuel  mixture.  Therefore,  the  author’s  calculation  data  are  true 
only  for  the  single  case  where  w,  a,  ®  ,  and  x  are  constant;  changes 
in  these  or  other  parameters  would  also  change  the  dependence  6^  «  fC^cs)* 

In  Yu.  A.  Shcherbina’s  work  [11]^  the  calculation  of  burnout 
behind  one  or  several  Ignition  sources  is  presented.  On  the  basis  of 
t-hat  calculation,  it  is  possible  to  construct  temperature  fields  at 
any  distance  from  ignition  sources.  It  is  assumed  that  the  flame 
front  is  oscillating  irregularly  in  relation  to  some  average  position 
a,  with  some  mean  square  deviation  d  from  that  average  position.  The 
value  of. ^  can  be  found  by  means  of  Taylor  equations  for  given  i  ,  s  , 
and  X. 


At  present,  the  value  of  u^  cannot  be  determined  analytically, 
and,  therefore,  the  author  has  found  an  empirical  dependence,  u»u  « 
f{ug,  w,  <y).  When  u^  is  known,  it  is  possible  to  find  the  value  a  **  f(x) 
and  consequently  to  calculate  the  temperature  profile  and  the  complete¬ 
ness  of  combustion  in  a  given  cross  section  of' the  combustor.  The 
following  observations  may  be  made  regarding  [Shcherbina’s]  work  [11]. 

As  a  matter  of  fact,  the  author  does  not  in  any  way  take  into 
account  in  her  calculation  the  effect  of  the  combustor  walls  on  the 
development  of  combustion  and  the  formation  of  the  flame.  Therefore, 
the  calculation  is  true  only  for  some  particular  cases  for  which  the 
dependence  a  ®  f(x)  was  found. 

Satisfactory  agreement  of  calculation  results  with  the  data 
from  experiments  by  other  authors  was  obtained.  However,  all  such 
experiments  apply  either  to  conditions  in  half-open  combustors  or  to 
c ompai’atively  short  combustors,  where  the  wall  effect  is  insignificant. 
Such  agreement  is  not  surprising,  since  the  author  employs  the  empiri- 
cal  dependence,  «  ^‘(uh,  w,  dj,  derived  from  these  same  experiments. 

Ac  tually,  she  considered  the  position  of  the  flame  in  the  combustor, 

=  f(x),  as  being  given,  when  in  fact  It  must  be  found  as  a  result 
of  the  calculation. 
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2.  Preliminary  Observations 


Before  turning  to  the  description  of  the  proposed  method  of 
calculating  the  combustor,  we  will  discuss  in  more  detail  some  char¬ 
s'  teristic  features  without  which  no  calculation  of  any  combustor  could 
be  caxTied  out.  Two  such  basic  features  are  the  velocity  uin^  of  tur~ 
tulent-flame  propagation  and  the  width  6^  of  the  combustion  zone  of 
a  .urbulent  fl^e  (or  the  value  i.e.,  combustion  time,  analogous 

in  meaning  to  ^t)* 


va 

to 

g' 

b-- 


i  C 


'  '  i 

6  X 
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At  present,  it  is  not  possible  to  calculate  theoretically  the 
lue  of  uip.  In  a  number  of  investigations,  attempts  have  been  made 
determine  experimentally  the  value  of  Ugi  in  an  open  flow  of  homo- 
"cous  fuel  mixtures.  The  question  of  how  to  determine  Urp  is  a  de- 
atle  matter  (see  [2],  [j,],  [9],  and  [10]),  an.d  since  the  results 
experimental  determination  of  u«p  depend  on  the  method  chosen  they 
considerably  with  different  authors.  In  cases  where  the  de- 
idtnces  u/p  •«  f(u*,  uh  p),  as  found  by  various  authors,  are  more  or 
. ?  slrail^  [i],  the  absolute  values  of  u^  differ  considerably  (thus 
(i  to  6)u').  A  change  in  the  absolute  value. of  affects  con- 
d-rably  the  results  of  calculation  of  and  r  *  fCx)-  Moreover, 
one  has  ever  mea  ired  u*r  in  any  combustor.  Therefore,  in  order  to 
this  problem,  we  made  an  attempt  to  determine  the  values  of 
utilizing  the  results  of  the  processing  of  experimental  investi- 


of  burnout  in  the  combustor. 


a  I 


The  matter  of  determining  the  ^  or  qp  is  still  more  contro- 


In  Talan^'''v'z  work  [8J,  there  is  a  derivation  of  the  theoret¬ 
ic  pendonct.  >  w,  w',  and  a  comparison  between  this 

niencc  [  Ic ]  and  experimental  results.  The  derivation  of  these 
:r^Lical  dcpcndenccz,  as  carried  out  by  the  author,  is  far  from 
,  .j  jel  f-evid'-iit ,  and  the  experimental  determination  of  gives 
to  a  number  of  objections. 

The  author  calculates  the  combustion  time  in  this  way.  He  con- 
;  the  combustion  of  an  individual  nucleus  of  fresh  mixture,  cut 
i  I /  turbulent  pulsation  from  the  initial  flame  surface  and  surrounded 
mbustion  products.  This  assumption  is  not  obvious.  Further, 

.  1  t.  tov  assume.-  that  such  a  nucleus  of  fresh  mixture  starts  burning 
■e  surface  with  a  speed  at  the  initial  moment  of 


—  “o  +  ~ , 
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where  Kn*  is  a  coefficient  In  which  the  effectiveness  of  pulsation- 
speed  action  is  taken  into  accotint  (the  value  of  this  coefficient  is 
not  determined),  w’q  and  1q  are  the  initial  values  of  the  pulsation 
speed  and  of  the  dimension  of  the  source,  respectively,  and  1  is  the 
running  dimension  of  the  nucleus,  which  burns  with  the  speed  Ujg  *  ug. 

Also,  no  allowance  is  made  for  the  increase  in  normal  velocity 
due  to  the  small  dimensions  of  the  source,  i.e.,  the  sharp  curvature 
of  its  surface. 

The  author  writes  the  equation  of  combustion  in  the  form 

dl’"U^X 


and  after  integration  he  arrives  at  the  final  form  of  the  equation: 


X  Bss-J^  In 


it  is  assumed  that  Kjj^  *=  1. 


Pig.  2.  Diagram  of  right  circular  cone  of  flame 
1  -  combustion  zone. 


The  entire  derivation  is  based  on  many  a  priori  assumptions, 
and  therefore  the  results  cannot  be  considered  as  proved. 

In  order  to  determine  ‘for>  the  width  of  combustion  zone 
was  measured  along  the  axis  of  the  right  circular  cone  of  the  flame. 
In  a  number  of  works  (e.g.,  [4]),  it  is  shown  that  the  value  of 
depends  on  the  distance  5c  to  the  ignition  source;  this  is  in  no  way 
accounted  for  by  Talantov.  Such  a  result  can  be  determined  from  an 
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clcsicntary  analysis  (seo  Pig.  2),  thun; 
by  definition  ,  M 


It.  is  possible  to  consider  approx liaatcly 


....ere  X  is  the  diatar.ee  froa  tr.e  fiocic-conc  apex  to  tr.c  base  of  t^.c 
tuiTicr,  d  is  \uo  clanetcr  of  the  burner,  and  Oij>  is  the  radial  Kldt.. 
of  tJx;  conbustion  zone.  Pros  this 


.  D  «  const.  ■  u* 


for  saall  x; 
for  large  x. 


(•> 

(••> 


Since  u<j.  doe.-;  not  depend  on  x  the  value  oc  in  both  erteo 
[{*)  and  (*♦).  d<'.  ends  cn  the  distance  to  the  ignition  sorcrce.  This 
;.clnt  is  not  considered  at  all  in  Talantov’s  work,  not  to  laention  the 
;  .  jt  that,  ger.eraily,  the  characteristic  of  Tqp  ia  quite  relative. 

In  ivallt: ,  the  cosbuctlon  tine  In  a  turbulent  flaac  Is 
-..teixalned  in  t.he  s.'jse  vay  cc  that  In  a  lanlnar  floroe: 


/ 


Acco.xllng  to  the  taodem  vlev.c  I’cgardlng  turbulent  cocbustlci!, 
otated  for  cxanplc  in  ll4j,  [15j»  and  otlwr  works  bas6d  on  Shchelkin'o 
,  ;-yalc.-l  model,  t^.c  width  of  t.ne  conbuotlon  zone  of  a  turbulent  flare 
Ic  defined  by  the  depth  of  curvaturH?o  in  the  flame  front  under  the  ac¬ 
tion  of  turbulent  pulsation.  Tnc  cu.-^atures  and  also  the  pulaatlcn 
o'’  tne  flame  fi'crit  ai'e  detemined  by  the  statla^cal  value  of  the  mean 
atatistlcal  displacement  of  the  flame  clement  vx*,  according 

to  what  was  said  above,  results  In  the  equality  6iji  -  vY*. 

If  no  allowance  is  made  for  the  effect  of  nomal  velocity 
ur.J  autoturbullzatlon  (as  a  xoault  of  smoothing  effect,  noixoal 
velocity  reduces  the  quantity  (Y*)''*;  this  turns  out  to  be  Justified 
in  a  number  of  cases),  the  quantity  ,  con  be  calculated 

ti.coretlcally: 

(ir’)  '  ^  IX  for  X  t'  . 

(h’)  for  »  > 
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Thus,  to  detemim  the  value  statistical  process¬ 

ing  of  a  great  number  of  instantaneous  positions  of  the  flams  front 
is  required. 

In  practice,  however,  the  detei^natlon  of  ttvs  width  of  the 
c<^ustion  zone  is  usually  carried  out  by  measuring  flame  temperature 
with  thermocouples  or  by  other  inertial  measursjsent  methods. 

The  data  obtained  by  means  of  ionization- type  pickups  along 
the  width  of  the  combustion  zone  may  be  of  definite  interest  because 
such  a  measures^nt  method  is  quite  sensitive  and  isaioss  it  possible  to 
register  individual  flame  oscillations. 

In  [4],  the  wiath  of  the  cc^ustlon  zone  was  determined  by 
ionization- type  picioips  rapidly  moving  across  the  flac^  (Fig.  3)  in 
an  open  flow  and  In  a  cylindrical  coaabustor,  speed  with  which 

the  flaise  was  crossed  amounted  to  abmit  0.5  m/sec,  while  the  velocity 
of  the  incident  flow  was  ^  %  m/&ec.  Actually,  not  aa,  instant neous, 
but  an  average- time  picture  of  the  distribution  of  ionization  current 
in  the  cross  section  of  the  flajE®  was  ta^n  in  this  way. 


Pig.  Determination  of  flame  Ijoundaries 
by  ionization  pickups 

1  -  homogeneous  gaaoline-alr  mixture;  2  - 
mixture- Igniting  source;  :  -  direction  of 
motion  of  an  ionization  pickup  across  the 
flas^,  4  -  typical  pictures  of  the  distri¬ 
bution  of  Ionization  current,  5  -  products 
of  complete  ccsabustion  of  the  mixture 

It  is  known  that  in  the  region  of  chemical  transformations, 
t.g.,  of  the  iaalnar  flame,  an  increased  ionization  car*  be  observed, 
it  it  just  this  phenomenon  which  makes  it  possible  to  sepai^te  the 
region  of  chemical  transformations  frcsa  bhat  of  the  products  of  cc®- 
;  letf*  combustion  in  a  turbulent  flame. 
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The  width  6ij^  of  the  coKbustion  zone  (along  the  radius  of  flame), 
detemined  in  [4]  on  the  basis  of  relations  from  the  probability  theory, 
should  (allowing  -1^  measurement  error)  amount  to 


(IG> 

In  order  to  verify  the  correctness  of  such  a  relation,  the  re- 
suits  of  experiments  from  (4j  were  processed  in  the  fom  of  the  re¬ 
lation 


K'  /(v).  where 


On  the  other  hand,  relations  (l4)  and  (15)  are  known. 

As  seen  from  Fig.  4,  the  quantity  6ip/4,6  found  and  the  values 

of  calculated  from  equations  (l4)  and  (15)  practically  coincide. 

The  velocity  of  the  incident  flow,  the  fuel-mixture  composition,  and 
therefore  u  j,,  did  not  affect  the  value  of  Om.  Also,  the  value  did 
not  change  as  a  result  of  the  specific  conditions  under  which  it  was 
determined  —  in  an  open  flow  or  in  a  cylindrical  combustor.  The  re¬ 
sults  obtained  have  also  been  confirmed  by  the  data  in  [11]. 


Fig.  4.  The  dependence  of  the  value  of  the  mean 
square  displacement  of  the  flame  on  the  relative 
distance  to  the  ignition  source 

— — —  -  calculated  values  of  the  mean  square  dis¬ 
placement  0  -  experimental  values  of  F*!*  found  in 
open- flow  installations  (W  ««  35  to  50  m/sec  and 
a  e  0.5  to  1.3);  e  -  values  of  “By  in  cylindrical  com¬ 
bustor  145  ram  in  diameter  (w  w>  35  to  75  iVsec  and 

a  »  0.7  to  1.3) 
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3. 


Proposed  Calculation  Method 


The  foregoing  discussion  indicates  that,  in  principle,  the 
calculation  of  a  consbustor  is  feasible.  As  was  done  in  a  nu^er  of 
the  works  discussed  above,  the  equations  (l),  (4),  (6),  and  (8)  can 
be  used  as  a  basis.  Also,  0.  Tzyan’s  condition  l9]  can  be  used,  as 
was  done  by  Talantov,  to  close  the  system  of  equations.  Finally,  the 
last  equation  needed  for  the  calculation  is  the  relation 


S-(n*=/(x) 

(see  the  equations  (l4)  and  (15) )• 

To  calculate  the  position  of  the  boundaries  of  the  flame  and 
of  the  burnout,  one  should  know  the  values  of  changes  in  velocity 
and  temperature  in  the  cross  section  of  the  flame  rather  than  the 
values  Cy  ^3  the  processing  of  experimental  results  shows, 

the  temperature  profile  along  the  cross  section  of  the  flame  is  quite 
complicated  (Fig.  5).  In  [11]  it  is  shown  that  the  temperatui^e  pro¬ 
file  can  be  determined  by  means  of  the  relationship 


whore 


'  WJI  •  •  /  0  i 


9 


®(t)  is  a  Gaussian  Integral  obtained  from  the  table. 

However,  such  a  law  of  temperature  change  makes  it  impossible 
to  integrate  equation  (6);  (a  nvimerlcal  integration  results  in  ex¬ 
cessively  cumbrous  equations). 

Analysis  of  the  existing  experimental  material  makes  it  possi¬ 
ble  to  find  a  more  simplified  solution.  The  law  of  temperature  change 
is  Introduced  in  a  linear  formj  then,  the  entire  change  in  temperature 
from  Tq  to  Tmax  takes  place  in  the  zone  at  the  depth 


The  coefficient  "3"  la  found  as  a  result  of  processing  many  ex¬ 
periments  (see  Fig.  4).  In  other  works,  e.g.  [6],  the  law  of  tempera¬ 
ture  and  velocity  changes  was  also  assumed  to  be  linear,  but  at  the 
depth  of  *•  6Gf). 
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Fig.  5.  Experimental  and  theoretical  dlmenBlonless 
temperature  profiles  behind  a  single  stabilizer. 

(A  linear  law  of  temperature  change  Is  assumed  In 
this  work:  solid  lines.) 


change  of  velocity  with  ia  also  assueied  to  linoas*} 


After  all  sii^plificationfij  the  system  of  equations 
tenaining  %t»  f(x)  appears  in  the  fora* 

*• 

(m 

f  1  —  t  J  , 

■■-)/ 
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!' 

(22) 

-  1 /jeL-i, 

“  %K  ^  y  i/.»  A  ’ 

(23) 

Where 

'f  -»ll  \  '• 

)‘ 

for  4t- 


4  t 


It  is  necessary  to  find:  *^0?  V 

fhe  Icsiown  values  arej  ic,  u^,  and  ug* 


There  are  eight  unknewn  values  and  seven  equations;  however, 
one  of  the  unknowns  is  an  argmaent  (e.g.^T^e  or^x).  Therefore, 
there  is  a  sufficient  number  of  equations  to  find  the  solution. 
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Fig.  6.  Pos.^tions  of  flame  boundaries  —  according 
to  Talantov’e  calculation  (l)  and  according  to  the 
calculation  by  the  proposed  method  (2).  [3  -  calcu¬ 

lation  results  by  the  proposed  method  with  * 

2.26  m/sec j  (the  curves  ”2”  were  found  for  w  « 

50  m/sec,  ujp  «  10  m/sec,  and  »  288*  abs.)3 


These  equations  make  it  possible  to  carry  out  the  calculation 
in  stages,  from  one  cross  section  to  the  next,  by  the  method  of  suc¬ 
cessive  approximations. 


To  compare  this  calculation  method  with  that  of  Talantov,  the 
example  given  in  [9]  was  calculated  for  the  following  initial  condi¬ 
tions:  w  •  50  m/sec,  uy  ■»  10  Vsec,  and  Tg  «*  288®  abs.  It  should 
be  noted  that  the  value  chosen  for  ut  in  the  example  is  considerably 
elevated.  Under  ordinary  conditions  at  s  »  0.05>  the  values  of  u^;* 
are,  as  will  become  evident,  noticeably  smaller  than  those  given. 

The  results  of  such  a  calculation  are  shown  in  Fig,  6.  It  can  be 
seen  that,  as  it  was  to  be  expected,  the  initial  positions  of  the 
flame  front  coincide  quite  closely  for  both  calculation  methods; 
however,  the  final  position  of  the  front  differs  considerably  fr<mi 
the  data  in  Talantov 's  calculation.  For  a  more  cozrect  comparison 
of  both  methods,  the  fact  that  the  values  of  uj  used  by  Talantov 
differ  considerably  from  the  real  one  should  be  taken  into  account. 

This  is  due  to  the  fact  that  um  was  determined  on  the  basis  of  the 
front  boundary  of  the  flame  ana  that  such  a  deteiislnation  is  incorrect 
(see  [2]).  In  reality,  as  follows  from  the  cflculation  based  on  our 
experln^ntal  data  (which  should  give  a  true  amount  for  um),  uj  »  2,26  Vsec 
for  e  *  0,05  and  a  »  1.2.  Under  the  same  conditions  ana  even  for  sosewhat 
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eBftller  of  Ug,  Talantov’s  •quatioa  [7]  gi^«8 

Uf  «  5.3(u’ )®‘ ’(^h)®**  ^  T.7  B/s«ee 

yfhm  tha  valuo  Mm  m  2,26  a/^^ee  is  ussd  and  th«  position  of 
fXsa®  boundaries  is  ealeulated  by  th«  proposed  aothod  and  than  co»« 
pared  with  the  ealculatlon  according  to  fiiXantov^  a  consldarsble 
discrepancy  in  the  results  is  seen  (see  Fig.  6). 

In  order  to  calculate  the  total  length  of  the  cwibustor,  it  is 
necessary  to  find  a  aethod  for  deterainlng  the  coordinates  of  nn  in 
the  combustor  cross  sections  where  «  0. 

In  this  case^  it  is  recosnended  that  the  following  self-evident 
relation  be  used: 


where  I«3.r 
rent  and  mq 
buation  son 


L 


is  the  width  of  the  ct^buation  zme  along  the  line  of  cur- 
)  is  the  average  flow  velocity  in  the  center  of  the  coa- 

^3.rl 


<p  2 


In  many  cases,  it  is  possible  to  assuae  approxiaately  u^p  ^  Ug; 
then,  the  final  e^qpression  is 


In  a  number  of  approxiaate  calculations,  when  the  shape  of 
the  curve  tjq  ■  f(x)  can  be  disregarded,  the  calculation  of  the  position 
of  flame  boundaries  [^c  «  f(x)}  becomes  substantially  simpler.  In 
such  a  case  the  total  length  of  the  cold  portion  of  the  combustor  (to 
the  point  where  the  flame  touches  the  combustor  walls,  1.®.,  at  i^c  * 
remains  practically  constant.  In  this  case,  the  entire  calculation 
can  be  carried  out  (byj  assuming  the  flame  front  to  be  infinitely 
thin,  l.e,,  by  finding  the  curve  nc  ■  f(x)*  *«d  then  finding  the 
value  6i«  •»  3®  for  each  running  value  of  x,  B:ie  calculation  Is  carried 
out  as  follows, 

Given:  1^,  k,  Xg,  t 

To  be  found:  xg,  uq  #  un  ,  nc 
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Of  ite  aS  fei>e  <iata  on  absolute  values 

values^  usiSTtL  ^5  toterestlag  to  calculate  the  u™ 

i Triigf  f  rn^i^co^ 

sssS^  al-ia^  ^"2S«  2.;rss'.sa.as“- 


«T 


r 


M 

[(t) 

1^1  -  a*,  -  (-^1  (1  ~  «ttj 

1’ 

a- 

(t)l 

(ncs+nei)’f4jc>  +  {%i~’)ci 

t)*l 

T„  =  ■  —  Z-iP  -f  Y ^p^Qi 


(24) 

(25> 

(26> 

m 

(28> 

(29) 

(30) 


of  neseScVlnXn^SS?  expenl»ental  results 

in  dlaitf?er*'Lfwf|f  nrhJtySSf  m''?n  ^5°  » 

and  for  stabilization.  StatS  preswel  Sfre  n?oJ»a 

bustor  (see  fig.  7).  ^  yfere  picked  up  along  the  jom- 
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Fig,  7.  Diagram  of  ooaaustor  measurement  system 
and  a2?rangement.  1  -  Position  of  turbulizatioa 
grating;  2  -  flame  staljilizer;  3  -  transparent 
windows;  4  -  flame;  5  -  statio  pressure  piekap 
points 


The  following  conditions  were  made  variable;  the  velocity  of  in¬ 
cident  flow  and  mixture  composition  a.  ©le  experimental  results  were 
found  in  the  form  of  the  function  z  »  f(z)  for  various  Initial  con¬ 
ditions. 


Since  the  value 


in  a  given  cross  section  of  combustor  is  Saiown  from  experiments,  it 
is  possible  to  determine  the  running  values  of  ^u,  Um,  et,c.  The 
values  of  dx  and  ^  are  known  (in  this  work,  it  was  assuisned  that 
s  »  0.05;  for  these  conditions,  the  value  of  o  is  detersiined  frm 
Pig.  4). 

The  values  of  u^i  were  determined  for  individual  portions*  dx 
of  the  combustor,  after  which  the  mean  value  was  calculated.  The  re¬ 
sults  of  the  calculations  are  shown  in  Figs,  8  and  9. 

The  values  found  for  u^  were  of  the  order  0.5  to  1,5  u‘;  this 
agrees  well  with  modem  views  about  turbulent  c^bustlon  (see  13.43 
and  [15]).  The  relationships  found,  u^  s  (u')®*  ®®  s^nd  Ug»  s  ug 

(in  the  range  a  ^  0.9),  are  in  satisfactory  agreement  with  known  ex¬ 
perimental  resiats,  which  give  «  f(u*,  U|^  ([153)#  taking  into 
account  calculation  errors. 

The  fact  that  the  effect  of  u*  (through  the  influence  of  w)  and 
uh  on  the  value  of  u^i  is  somewhat  stronger  than  that  in  [153  o®n  be 
explained  by  the  effect  of  the  ignition  source.  In  the  case  of  a  small 
ignition  source,  the  cos^leteness  of  ooxobustion  close  to  the  ignition 
source  in  the  axis  of  flow,  as  is  known,  may  be  considerably  less  than 
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Pig.  8.  Relative  change  in  the  velocity  of  tur¬ 
bulent  combustion  in  relation  to  that  of  normal 
combustion 
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Dependence  . . . 

^T’.T0K 


♦  2  -  dependence 

^H.TSK 


Um  Uf*  ^2 

-  «  —2 -  X.  - -  ;  3  -  dependence  -  cal 

^T.Ten  ^K.TeK  '^2t8k 

culation  results  from  experimental  data  for  Wj,  « 
50  m/sec,  a  »  1  1.6  [TEK  «  running] 


Pig.  9.  Relative  change  in  the  velocity  of  turbu¬ 
lent  combustion  with  increasing  flow  velocity 
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the  usual  0.95  to  0.98.  In  this  case,  the  Ignition  of  the  fuel  mix¬ 
ture  does  not  start  close  to  the  stabilizer  edges  but  further  down 
the  flow;  this  is  not  taken  into  account  in  the  calculation. 

It  should  be  noted  that  since  under  the  majority  of  operating 
conditions  the  flame  front  does  not  touch  the  coaabustor  walls  $^ft<l,0), 
corrections  for  expansion  of  combustion  products  were  Introduced  when 
computing  u^  during  the  calculation  of  u^  w  f(uH)>  was  done  in 
Yu.  A,  Shcherbina’s  work  [11]. 

Thus,  finally,  in  the  first  approximation  it  is  possible  to  take 


// 


T 


Where  Ut©  *=  O.092,  Wq  «  30  m/sec,  **  8l  cm/sec,  Tjaax© 
and  Sg  «  0.05. 


2400®  abs.. 


Fig.  10.  Position  of  flame  boundaries 

1  -  from  calculation;  2(o)  -  experi¬ 
mental  points  determined  by  the  ioniza¬ 
tion  pickup  method;  3  -  boundaries 
found  by  direct  photography 
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Pig.  11.  Position  of  flame  boundaries 

1  -  from  calculation;  2(o)  -  experiment¬ 
al  points  determined  by  the  ionization 
pickup  method. 


Pig.  12.  Position  of  flame  boundaries 

1  -  from  calculation;  2(o)  -  experiment¬ 
al  points  detemined  by  the  ionization 
pickup  method. 


In  the  case  of  the  flame  touching  combustion  \7all8, 

/’‘■•"in 
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In  addition*  a  comi^ari con  v/a^  made  of  tht  restiitc  of  the  de- 
*:<  mination  .lamt  boundaries  in  a  cylindrical  e<^buiitor  by  the 
calculation  mt^.aod  witn  those  detersiined  by  the  Ionisation  pickup 
method  [rj*  and  alco  b'j  photograph^"  through  a  windot^?  In  the  com- 
biistcrn  As  c,ecr.  from  Figs.  It.  to  12,  good  agreement  with  calculation 
iata  wa.  found. 


Cc-nciusion 


A  calculation  Hitthod  for  tht  .  liaplest  ramJet-type  combustor 
i.  :rc..<.ntcd.  The  method  con.^i..;ts  in  closing  the  system  of  the  basic 
ga. -dynamic  equations  (of  fuel-mixture  flow)  within  a  cylindrical 
tube  with  a  point- Igidtion  :ource,  by  introducing  an  auxiliary  con¬ 
dition  which  defines  the  dependence  of  the  width  of  the  combustion 
June  on  the  di3tai»co  to  the  ignition  source.  The  results  of  the 
calculation  wore  compared  with  some  experimental  data  and  were  found 
to  be  In  good  agreement  with  them. 
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